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Abstract 

The relationship between Elhptic Ruijsenaars-Schneider (RS) and Calogero-Moser 
(CM) models with Sklyanin algebra is presented. Lax pair representations of the 
Elliptic RS and CM are reviewed. For n = 2 case, the eigenvalue and eigenfunction for 
Lame equation are found by using the result of the Bethe ansatz method. 
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1 Introduction 



A general description of classical completely integrable models of n one-dimensional particles 
with two-body interactions V{qi — qj) was given in Ref. 1). To each simple Lie algebra and 
choice of interaction, one can associate a classically completely integrable system^^"^' such 
as a rational, hyperbolic, trigonometric or elliptic CM model. 

The Lax pair representation (Lax representation) of a system is a direct method of 
showing its integrability and the complete set of integrals of motion can also be constructed 
easily. The Lax representation and its corresponding r-matrix for rational, hyperbolic and 
trigonometric A^-i CM models was constructed by Avan el al.^^ The Lax representation for 
the elliptic CM models was constructed by Krichever^^ and the corresponding r-matrix was 
given by Sklyanin^^ and Braden et al.^^ There exists a specific feature in that the r-matrices 
of the Lax representations for these models turn out to be dynamical (i.e., they depend on 
the dynamical variables) and satisfy dynamical Yang-Baxter equations. ^^'^^'^^'^^ 

For the dynamical r-matrix, the fundamental Poisson algebra of the Lax operator, whose 
structural constants are given by a dynamical r-matrix, is generally no longer closed. The 
quantization problem and its geometrical interpretation are also difficult. Considering all of 
these, a non-dynamical r-matrix is found for these systems. ^^-"'^^^ The trigonometric limit of 
these results can be found in Ref. 42). We know the Lax representation for a completely 
integrable model is not unique. The different Lax representations of an integrable system are 
conjugate to each other (for the field system they are related by a 'gauge' transformation). 
The corresponding r-matrices are related by a 'gauge' transformation which is the classical 
dynamical twisting relation^^^ between those r-matrices. 

The RS model is a relativistic generalization of a CM model. It describes a completely 
integrable system of n one- dimensional interacting relativistic particles. It can be related 
to the dynamics of solitons in some integrable relativistic field theory. ^^■^^^'^^'■'^^^ Its discrete- 
time version has been connected with the Bethe ansatz equation of the solvable statistical 
model. Recent developments have shown that it can be obtained by a Hamiltonian reduc- 
tion of the cotangent bundle of some Lie group, and can be considered as the gauged WZW 
theory. The Lax representation and its corresponding r-matrix for rational, hyperbolic and 
trigonometric A^-i type RS models was constructed by Avan et al.^^ The Lax representa- 
tion for the elliptic RS models was constructed by Ruijsenaars,^^) and the corresponding 
r-matrix was given by Nijhoff^^^ and Suris.^^^ The main difference between the r-matrices of 
the relativistic (RS) and non-relativistic (CM) models is that the latter is given in terms of 
a linear Poisson-Lie bracket, whereas the former (RS model) is given in terms of a quadratic 
Poisson-Lie bracket. In contrast with the dynamical Yang-Baxter equation of the r-matrix 
for the CM model. the generalized Yang-Baxter relation for the quadratic Poisson-Lie 
bracket (RS model) with a dynamical r-matrix is still an open problem.^°^ Moreover, the 
Poisson bracket of the Lax operator is no longer closed, and consequently the quantum ver- 
sion of the classical L-operator, has not been constructed. However, as for the CM model, 
a different Lax representation which is conjugated to the original one can be found. The 
corresponding r-matrix changes by a 'gauge' transformation. The resulting r-matrix may 
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be non-dynamical. Such a transformation may be called the classical dynamical twisting of 
the associated linear Poisson-Lie bracket. Due to the quadratic Poisson-Lie bracket of the 
RS model, there exist dynamical twisting relations between the r-matrices of Lax operators 
related by gauge transformations. Such dynamical twisting is the semi-classical limit of the 
quantum dynamical twisting of the i?-matrix in Ref. 12). For recent progress in the study 
of CM models, see, for example, Refs. 21)-24). 

The paper is organized as follows: In Sec. 2, we present some general formulae for 
dynamical systems. In Sec. 3, we review some results for the elliptic RS and CM models. 
The non-dynamical r-matrices for the integrable elliptic systems are then presented. Their 
quantization conditions correspond to the quantum Yang-Baxter relation, and the i?-matrix 
is simply the Z„-symmetric Belavin model.^^^ In Sec. 4, we will present the relationship 
between the Sklyanin algebra^^'^^^ and the integrable systems. In Sec. 5, we will obtain the 
eigenvalue and eigenfunction for the Lame equation. The Lame operator is equivalent to the 
Hamiltonian of the elliptic CM model. Section 6 has some brief summary. 



A Lax pair (L, M) consists of two functions on the phase space of the system with values in 
some Lie algebra g, such that the evolution equations may be written in the following form 



where [,] denotes the bracket in the Lie algebra g. If we have formulated the Lax pair 
relation, the conserved quantities (integrals of motion) can be constructed easily. It follows 
that the adjoint-invariant quantities trL'(/ = l,...,n) are the integrals of the motion. In 
order to implement the Liouville theorem onto this set of possible action variables we need 
them to be Poisson-commuting. As shown in Ref. 25), the commutativity of the integrals 
trL' follows if the Lax operator can be deduced from the fundamental Poisson bracket 



2 The dynamical twisting of classical r-matrix 



dL 



[L,M], 




{L,{u),L^{v)} 



(2) 



or quadratic form^^^ 



{L,{u),L2{v)} 



Li{u)L2{v)r-^2{'^,v) - r^-i{v,u)Li{u)L2{v) 
+Li{u)s'^{u,v)L2{v) — L2{v)s~{u,v)Li{u), 



(3) 



where we use the notation 



L 



1 — 



1/(8)1, 



L2 = 1 ® L, 



021 — PdYiP, 



(4) 
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and P is the permutation operator such that Px ® y = y ® x. 

For the above relations to define a consistent Poisson bracket, one should impose some 
constraints on the r-matrices. The skew-symmetry of the Poisson bracket requires that 



ri2{u, v) - stiiu, v) = r^^iu, v) - Si2{u, v). (6) 

For the numerical r-matrices r^(M, v), sf2{u,v), some constraint conditions (sufficient condi- 
tions) are imposed on the r-matrix in order to make it satisfy the Jacobi identity. However, 
generally speaking, the r-matrices r^(w, v), v) depend on dynamical variables, and the 
Jacobi identity which implies an algebraic constraint for the r-matrices takes a very compli- 
cated form 

[Li, [ri2, ris] + [ri2, ras] + [r32, ^13] + {-^2, ^13} - {-^3^12}] + cyc.term = 0. (7) 

It should be remarked that for a given integrable system, we can choose different Lax 
formulations. The r-matrices corresponding to different Lax formulations are generally dif- 
ferent. So, in some cases, we can transform a dynamical r-matrix into a non-dynamical 
r-matrix. ^"^'^^^ The different Lax representations of a system are conjugate to each other: 
if (L, M) is another Lax pair of the same dynamical system conjugate to with the old one 
(L, M), it means that 

L{u) = g{u)L{u)g ^{u), 

M{u) = g{u)M{u)g-\u)-{j^g{u))g-\u), (8) 

where g{u) & G whose Lie algebra is g. Then we have 

Proposition: The Lax pair (L, M). has the following r-matrix structure 

{Li{u),L2{v)} = [fu{u,v),Li{u)] - [f2iiv,u),L2{v)], (9) 

where 

ri2{u,v) = gi(u)g2{v)ri2(u,v)gi^{u)g2^{v) + g2(v){gi{u), L2{v)g^^{u)g2^(v) 

+ l[{gi{u),g2{v)}g^\u)g2\v),g2{v)L2{v)g2\v)]. (10) 
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For a given Lax pair (L, M) and the corresponding r-matrix, if there exists a g such that 



h 



12 



9i{u)g2{v)ri2{u,v)g^ ^{u)g2 ^{v) + g2{v){gi{u), L2{v)]g^ \u)g2 \v) 
+ ^[{9iii^),92{v)}gi\u)g2\v),g2{v)L2{v)g2\v)] 



(11) 



and 



dn,h 



dpjhi2 







(12) 



then a non-dynamical Lax representation of the system exists. 

By a straightforward calculation, we can also find that the twisted Lax pair (L, M) and 
the corresponding r-matrix fi2 satisfy 

[Li, [fi2, fia] [fi2, f23] + [?'32, na] + {1,2, fia} - {^3^12}] + cycl.term = 0. 
Similarly, for the quadratic form, the Lax pair (L, M) has the following r-matrix structure 



where 



{h{u)Mv)} 



Li{u)L2{v)fi2{u,v) - f^2{u,v)Li{u)L2{v) 
+Li{u)st2{u,v)L2{v) - L2{v)si2{u,v)Li{u), 



(13) 



r^2[u,v 
rt2{u,v 

st2(u,V 
Si2{u,V 

Ai2(u,v 



gi{u)g2{v)r^2{u,v)g^^{u)g2^{v) - A^2{u,v) + A2i{v,u), 

gi{u)g2{v)rt2{u,v)g^\u)g2\v) - Ai2{u,v) + A2i{v,u), 

gi{u)g2{v)st2{u,v)g^\u)g2^{v) - A2i{u,v) - A^^^{u,v), 

gi{u)g2{v)s^2{u,v)g^\u)g2^{v) - Ai2{u,v) - A2i{v,u), 

L^\v) A12 {u,v), (w,v) = Ai2{u,v)L2\v), 
g2{v){gi{u), L2{v)}g^^{u)g2^{v) 

+7:[{9iiu), g2{v)}gi\u)g^\v), g2{v)L2{v)g^\v)]. 



(14) 



There are still relations: 



-ri2iu,v) 



rt2{u,v) - st2{u,v) = r^2{u,v) - Si2{u,v). 



(15) 
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And also we have, for given Lax pair (L, M) and the corresponding r-matrices, if there 
exists a g such that: 



9i{u)g2{v)st2{u,v)g^\u)g^^{v) - A2i{v,u) - A'ij{u,v) = 0, 
gi{u)g2{v)s^2i'^,v)g^\u)g2^{v) - Ai2{u,v) - A^i^(w,m) = 0, 
hi2{u,v) = gi{u)g2{v)r^2i'<J','")9l^i'<J')92^{v) - Ai2{u,v) + A'^2i{v,u) 

- gi{u)g2{v)rt2{%v)gi\u)g2\v) - A^^{%v) + aS{v,u) (16) 

and 



dqihu = dpMi2 = 0, 



(17) 



then a non-dynamical Lax representation with Sklyanin Poisson-Lie bracket for the system 
exists. 



3 Lax pair for elliptic RS and CM models 

We first define some eUiptic functions: 



^9 



1 _ i 

2 n 
1 

2 



a 
b 



u, nr), a{u) — 9 

, r) = E~ -oo exp {i7i[{m + afr + 2{m + a){z + b)]) , 
9'^\u) = du(9^^){u)),a'{u) = d^{a{u)),i{u) = d^{\na{u)), 



(18) 



where r is a complex number with Im(r) > 0. 

The Ruijsenaars-Schneider model is a system of n one-dimensional relativistical particles 
interacting by a two-body potential. In terms of the canonical variables pi,qi {i = 1, ■ ■ ■ ,n) 
enjoying the canonical Poisson bracket 



{Pi,Pj} = 0, {qi, Qj} = 0, {qi,Pj} = Sij, 
the Hamiltonian of the system is expressed as^^^ 



(19) 



H = mc^ ^ cosh pj Y[ 
j=i \ k^j 



(20) 
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where qjk = qj — qu, fn denotes the particle mass, c the speed of hght, and 7 is the coupling 
constant. The above defined Hamiltonian is known to be completely integrable.^^^'^^^ As we 
mentioned above, the Lax representation (Lax operator of the classical L-operator) is one of 
the most effective ways to show that the system is integrable. One Lax representation for 
the elliptic RS model was first formulated by Ruijsenaars:^^^ 



Here, we use another Lax representation L^^^ 



Lj^ can be obtained from the standard Ruijsenaars' Lr{u) by using a Poisson map 



Qi — Pi — ^Pi+olnl|— -. (23) 

2 a{q,k - 7) 

Following the work of Nijhoff et al.,^*^) the fundamental Poisson bracket of tb Lax operator 
Lr{u) can be given in the following quadratic r-matrix form with dynamical r-matrices 



{Lr{u)i,Lr{v)2} = LR{u)iLR{v)2r\2{u,v) -rt2{u,v)LR{u)iLR{v)2 

+Lij(M)is+ l'"' v)Lr{v)2 - Lr{v)2S^2{'^-: v)Lr{u)i, (24) 

where 

r\2iu, v) = ai2{u, v) - Si2(m) + S2i{v), r^^i^^ ^) = «i2(m, + ^1*2 + ^12^ 

st2{U:V) = Si2{u)+Ut2: S'[2{u,v) = S2l{v) - U^2: (25) 

and 



u 



± 
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ai2{u,v) = r'l2{u:v)+J2^{u-v)eii<^eii + ^^{qij)eii<^en, (26) 

1=1 i^j 

r?2(ii, v) ^ -7^^7—3-4^^ ® eji, Si2(u) = Y^{LR{u)d^LR{u))eij eji. (27) 
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The following properties are satisfied: 



rti {v,u) = -rt2 {u,v), sti{v,u) = 3^2 {u,v), 

ri2 {u,v)- st2 {u, v) = {u, v) - s^a {u, v) . (28) 

Here we would like to reformulate the Lax formulation for the RS model. Define an n <S> 
n matrix A{u] q) as: 

A{u; q)] = A{u, gi, ■ ■ • , = e^'\u + nq, -Y.(Ik + -^)- (29) 

k=i ^ 

Here we should point out that A{u\ q)j corresponds to the intertwiner function of cpj'^ between 

the Z„-symmetric Belavin i?-matrix^^^'^^) and the A^li face model. ^'^^'^'^^ 
Define 



g{u) = A{u;q)A{q), A{q)] = hi{q)6], 
hj{q) =hj{qi,---,qN)^ — ^ZTT" 
We can construct the new Lax operator L{u) as 



L{u)^g{u)LR{u)g-\u). (31) 
More explicitly, it can be written as: 



n 1 

= E -T^M^ + nr, qykA-\u; g)Je^"= ,i, j ^ 1,2, ■ ■ ■ ,n. (32) 

k=i ^(7) 

It can be proved that the fundamental Poisson bracket of L{u) can be given in the quadratic 
Poisson-Lie form with a nondynamical r-matrix: 



{Li{u),L2{v)} = [ru{u-v),Li{u)L2{v)]. (33) 
Here the numerical r-matrix is the classical Z„-symmetric r-matrix.^^^It takes the form 



r^^(v)=\ ^ "i)eit-i){v)e(i-t)io) + {e^i-^)iv) W)) ^/ ^ + J - ^ + mod n ^ ^^^^ 
I otherwise 
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We know the Z„ symmetric r-matrix satisfies the nondynamical classical Yang-Baxter 
equation 

[^12(^^1 - -^2), ^13(^1 - V3)] + [ri2(vi - V2), ^23(^2 - ^^3)] 

+ [ri3{vi - ?^3), ^23(^^2 - V3)] = 0, (35) 

We also know that this r-matrix is antisymmetric and Zn symmetric: 



Antisymmetry : —r2i{—v) — ri2{v), 
Zn (8) Z„ Symmetry : ri2{v) — {a<^ a)ri2{v){a (S> a)~^, (36) 

where a = g, h, and g, h are n®n matrices defined as: 



hij - (^j+ij mod n' 9i3 - ^'kj- (37) 
For convenience, we also define another n<S>n matrix 

Ia = Iaua,=9'''h"\ (38) 

where ai,a2 G Zn and u = exp{27ii / n) . 

Next, we will consider the non-relativistic limit of the Lax operator L{u). First rescale 
the momenta {pi}, the coupling constant 7 and the Lax operator as follows:^^'' 

Pi:^-Pp[, n^:^Ps, L{u) -.^ a{^)L' {u). (39) 

77/ 

Here notation L' is introduced. The non-relativistic limit is then obtained by taking the 
limit /5 — > 0. We have the following asymptotic properties 

L'{u)) ^8]- (3{Y.{A{u; q)lA{u; q)% - sd,^A{u; q)l)A-\u; g)J}) + 0{P'). (40) 

k 

If we make the canonical transformation 



P'i^P'i-Z;^^^M{q), M(g) = n^(%)' (41) 
It' uqi -^j 

we finally obtain the Lax operator of the eUiptic An-i CM modeP°) 
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L'(uY. — (5* 

Lcm{u)) = -lim j-^U-^v\-i^MMi,) (42) 



Here we have 



{Lcm{u)i, Lcm{v)2} = [ri2{u - v), Lcm{u)i + Lcm{v)2\- (43) 

For the newly constructed Lax representation L{u), the quantization becomes no longer 
difficult. Define the Z„-symmetric Belavin's i?- matrix as: 

r e'W(o)a(n)a(V— ft) e(o)(^)e('-^)(t,+v3TR) „•,„•_/, l „ 
1 otherwise. 
We know this i?-matrix satisfies the quantum Yang-Baxter equation 

Rl2{ui - U2)Ri3{Ui - %)i?23(M2 " %) = ^23('«2 " U3)Ri3{Ui - U3)Ri2{Ui - U2). (45) 

The i?-matrix is symmetric in the sense that 

Ri2{u) = {a ® a)Ri2{u){a ® a)~^ , a = g,h. (46) 
By taking the special limit ^ ^ 0, we can obtain the classical symmetric r- matrix 

Ri2{u)\n-^o =1^1 + ^/^hri2{u) + o{h^). (47) 
Now let us study the quantum L-operator, using the usual canonical quantization procedure 

, — d 

Pj ^Pj ^ -^-'^^0^,^ Qj^Qj^ j = (48) 
The corresponding quantum L-operator can be formulated as: 

1 " 

^iV k=i 



= — -E^(« + ^7;5)r^-'(«;g)fe-^^^^. (49) 

It should be remarked that this quantum L-operator is simply the factorised difference 
representation for the elliptic L-operator. ^^^'^^^ The above defined quantum L-operator sat- 
isfies the quantum Yang-Baxter relation 

Ruiu - v)Li{u)L2{v) = L2{v)Li{u)Ri2{u - v). (50) 
The proof can be found in Refs. 31), 34), 33) and 35). 
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4 RS and CM models related with Sklyanin algebra 

We introduce here some notation for elliptic functions: 



aaiu) = 9 



1 I ai 

2 n 

1 I 02 

2 ~ n 



-in) aQ{u + ^f-ih) 

The above mentioned quantum it!-matrix can be rewritten as following up to a scale: 



(51) 



(52) 



as before a = ai, a2 and ai G Zn, i = 1,2. 

The quantum L-operator L obtained in the last section can be represented by the gen- 
erators of Sklyanin algebra Sq,: 



L{u) = J2Va{u)IaSa, 



(53) 



where 



2 ' 



(54) 



where 5 is a constant. 

The quantum Yang-Baxter relation (49) gives the defining relations of the Sklyanin 

algebra:^^)'^^ 



^ C^7i-«i + (ft-7i)(7i-^2)^^^^^^_^^^^(0)^^(2v/^n) _ 



(T^_a ( V^^)0-a+/3-7 ( V^^) (7^{^J-lK) (T/3_^ ( V^^) 

with Qij, 7j e i = 1, 2. 

We can give a realization of the generators of Sklyanin algebra as: 



(55) 



(56) 



Here we introduce the symbol Sja- for the elliptic function 



(57) 
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Next, we will consider the classical limit of the above defining relations. Letting h ^ 0, 
the quantum i?-matrix become the classical r-matrix, and we explicitly have the elements 
of the r-matrix in (34) presented in the last section, here we use another notation 

R(u) = 1 + y/^hr(u) + 0{h^). (58) 
The classical r-matrix is written as: 



r(ii) = X]^«(^i)4«)/a^ (59) 

a 

where 



Wq{u) = 0, 

aa(?i)ao(0) , . 

Wa[u) = a^O. 60 

aa{0)ao{u) 

In order to consider the classical limit of L, we first present the classical hmit of Va{u): 



Vo{u) 
VJu) 



naa{0)ao{u') nao{u') 



(61) 



(62) 



From the definition of the operator Sa, we easily have 



(63) 



In the limit ^ — > 0, the elliptic functions take the forms: 



SjQ — 



q. — 

>~>]o. — 



-ih) 



i + v^^eE 



(64) 



(-i)-..(o) n ^ 



1 + 



^^a(0) t^^aaiqjk) 



(65) 
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So, we have 




Mljk) dqj 



-n 



d_ 

dqj 



nao{u') Yfe/j- MQjk) 



+ 0{h% 



(66) 



n- 



(67) 



here we have Y^i y2k^i '^oiij^) _ g because is an odd function. 

We can finally expand the quantum L operator in the order of h when we take a limit 
^ — > 0. However, we first introduce some notation 



L{u) = ^ V^{u)SJ^ = 1 + ^hl{u) + 0(^2), 



(68) 



where l{u) is the classical / operator. We may represent l{u) in terms of generators of the 
"classical" Sklyanin algebra Sa- 



l{u) 



nao{u') 



7T - ^Va{u)Sa. 



(69) 



The function Va{u) is defined as: 



vo{u) 



1 

n 



1 (Ta(li)(7a(0) 
^^a(M) = ^-7 — , a / 0. 



(70) 
(71) 



From the above obtained results, we can realize the generators of the "classical" Sklyanin 
algebra in the following forms: 



So = 



Srv — 



^ d 



(^a{q. 



'jk) 



-ir<-a(o)En 



(72) 
(73) 
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On the other hand, here we can say we give a definition of the generators of the "classical" 
Sklyanin algebra. 

In the classical limit, the quantum Yang-Baxter relation becomes the following 



[h{u),l2{v)] = [ri2(^-^;),/i(^)+,/2(^;)] (74) 

Substitute the /-operator with the generators of the " classical" Sklyanin algebra, and through 
tedious calculation, we have 

[s^, s^] = (u;"^^^ - u;°^^^) (^^^ (75) 
On the other hand, we find that 1^ satisfy a similar relation 



[7„,7^] = (u;-i^^-a;-^^i)7„+^. (76) 

So, after rescaling Sa , we find {sa}- and {la} satisfy the same algebra. 

Finally we should point out that if we substitute by the corresponding canonical 
variable pk, the /-operator will become a T-operator, and the commutative bracket on the 
left-hand side of the above relation (74) will change to the standard Poisson-Lie bracket. 
Here we rewrite as: 



{T,{u),T2{v)} = [ru{u-v),T,{u) + T2{v)]. (77) 

5 CM model, Gaudin model, Lame equation and the 
Bethe ansatz 

For the difference factorized operator L, we can find some commuting families which are 
related to conserved operators. By using the fusion procedure, the commuting family take 
the form 

L'„ = tr[L(„)®---®L(„)P™], 

there are m L s above, P!" is the anti-symmetric projector. In the classical limit, we also 
have a similar commuting family 
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where a* e Z^, i — 1, • • • , m. Let uf — 0, so u — uq = — nh^, and after rescaling a/(M), 
we have 



QmK) = E • • • SaJr[Ia, ® 4^^™], (79) 

We will discuss a special case n = 2, explicitly we have 

V 0-0(^12) cro{qi2) oqi dq2 J 
sio = 2aio[U) 4 — - — r — — j , 

V cro{qi2) cro{qi2) oqi oq2 J 

sn = 2.n(0) - '-^(f - ^)) . (80) 

V cro(?i2) (7o(gi2) dqi dq2 J 

We will calculate the non-trivial conserved operator 
After some tedious calculations, we have 



«2 = -r^ + 2C^(TT^-TT^)+4 



C^o(?) 0-o(?) 5?2 55l9?2 

where q = qi — q2- This relation is just the same as that obtained by Hasegawa.^^^ 

Since 1 and ^ + ^ are also conserved quantities defined above, after some tedious 
calculations we have another conserved operator 

We can change it to a more familiar form. Let ^ = 2/3, and suppose ip and A are an 
eigenfunction and eigenvalue of the above Hamiltonian 



Hxp = K^. (84) 
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At the same time, we introduce a transformation of this eigenfunction -0 = ip[ao{q)]^, we 
thus have the following relations: 



This means: 

H'i^ = + + 1)^ In '^o(g)]^ = A^. (86) 

One finds that H' is simply the Hamiltonian of the CM model, see, for example, Refs. 3), 
4), 37) and 38). It is also connected to the Lame operator, see Ref. 39) and the references 
therein. 

Next, we will calculate the eigenfunction and eigenvalue of the above defined Lame oper- 
ator H. Here wc first review some of the results obtained by Felder and Varchenko.^^-* The 
difference operator L which is equivalent to 5*0, one generator of the Sklyanin algebra when 
n — 2, is given by 

ao{q + 2h(3) ao{q-2hf3) 
Lil^yq)^ r\ V'(?-2^) + t-. ^{q + 2h). (87) 

This difference operator is also called the g-deformed Lame operator. In the frame work of 
the quantum inverse scattering method, there is a result as foUows:^^) 
Let {ti, - ■ ■ , tm, c) be a solution of the Bethe ansatz equations: 

ao{f, - 2hff) ao{fj - - 2/i) ^ 

ao{U + 2h(3)}},ao{tj-U + 2h) ' ''"'^^^ ^^^^ 

such that ti ^ tj mod Z + tZ if i ^ j. Then 

P 

^{q)^e^llaoiq + tj), (89) 
j 

is a solution of the g-deformed Lame equation Lip — eip with eigenvalue 

^ ,n Mm A Mtj + (2/3 - 2)h) 

ao(2h(3)j-}, ao(tj + 2(3h) " ^''''^ 

By taking the special limit /i — > 0, the difference operator becomes: 
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L = 2 + 4^^ 



dq 



ao{q) dq ao{q) 



+ o{h'). 



(91) 



We find the term with order h"^ is exactly the Hamiltonian presented in relation (85) 



cT^(g) d 



dq"^ 



ao{q) dq ao{q) ' 



(92) 



Since we already know the eigenfunction of the difference operator L is ip{q) — e"^^ 11^=1 co(?+ 
tj) which docs not depend on h, we need only expand the eigenvalue of L in the order of h. 
We can obtain the eigenvalue of the Lame operator H 



e = 2-4h 



1\ (^o{tj) 



^o(^i) ( ^o(^i) 



(To{tj)_ 



i>j 



cro{ti)ao{tj) 



i=i 



+ 0(^1^). 



At the same time we take the limit ^ ^ for the Bethe ansatz equation, obtaining 



(93) 



c + p 



(94) 



Considering this Bethe ansatz equation, we can finally find the eigenvalue of the Lame 
operator A, 



A = (1 - 2/3) 



2^0 (0) 



(95) 



Here we have the results: 

Let (ti, • • • , tjn-i c) be a solution of the Bethe ansatz equations: 



(Jo(ii) j.^j-^^ (Jo(ii - tj) 



0, ^ = l,---,/3. 



(96) 



such that ti 7^ L mod Z + rZ if i 7^ j. Then 
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i^{q)^e'^'J{<J,{q + t^) (97) 

3 

is a solution of the equation 



dq 



ao{q)dq ao{q) 



ip{q) = Aip{q), 



(98) 



with eigenvalue 



A= (1-2/3) 



, ^2^0 (0) 

^ <^^(0)' 



(99) 



We can also obtain these results by directly using the algebraic Bethe ansatz method for the 
A^i^ Gaudin model,^'^) just like the algebraic Bethe ansatz for the XYZ Gaudin model.^^^ 



6 Summary 

We review some developments concerning the non-dynamical structure of the elliptic RS and 
CM models. We also give a solution to the Lame equation. The eigenfunction and eigenvalue 
for the Lame operator are found through the results of the Bethe ansatz. 

The results of the last sections are only for the n — 2 case. For general n, we can 
also obtain the eigenvalue and eigenfunction for the generalized Lame operator by using the 
algebraic Bethe ansatz method for the sl{n) elliptic Gaudin model. The conserved quantities 
also correspond to the Hamiltonian of the elliptic CM model. 
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